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ABSTRACT 



While the definitive detection of gamma-rays from known 
supernova remnants (SNRs) remains elusive, the collec- 
tion of unidentified EGRET sources that may be associ- 
ated with SNRs has motivated recent modelling of TeV 
emission from these sources. Current theoretical models 
use power-law shock-accelerated protons and electrons in 
their predictions of expected gamma-ray TeV fluxes from 
those unidentified EGRET sources with remnant associa- 
tions. In this paper, we explore a more detailed non-linear 
shock acceleration model, which generates non-thermal 
proton distributions and includes a self-consistent deter- 
mination of shock hydrodynamics. We obtain gamma-ray 
spectra for SNRs allowing for the cessation of acceleration 
to high energies that is due to the finite ages and sizes of 
remnants. Gamma-ray spectral cutoffs can be observed 
in the TeV range for reasonable remnant parameters, and 
deviations from power-law behaviour are found at all en- 
ergies ranging from 1 MeV up to the cutoff. Correlated 
observations by INTEGRAL, Whipple and other instru- 
ments may provide stringent constraints to our under- 
standing of supernova remnants. 
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1. INTRODUCTION 



Supernova remnants have long been invoked as a princi- 
pal source of galactic cosmic rays, created via the pro- 
cess of diffusive Fermi acceleration at their expanding 
shock fronts (e.g. Drury 1983, Lagage and Cesarsky 1983, 
Blandford and Eichler 1987). Such systems can also pro- 
vide gamma-ray emission via the interaction of the cosmic 
ray population with the ambient remnant environment; 
this concept was explored recently by Drury, Aharonian 
and Volk (1994). In their model, the gamma-ray lumi- 
nosity is spawned by collisions between the cosmic rays 
and nuclei from the ambient SNR environment. Inverse 
Compton scattering involving shock-accelerated e~ and 
the cosmic microwave background and also IR/optical 
emission (from dust /starlight) forms added components 
in the 7-ray SNR models of De Jager and Mastichiadis 



(1996) and Gaisser, Protheroe and Stanev (1996). 

Given the absence of definitive detections of gamma-rays 
from known supernova remnants, the motivation for mod- 
elling these "hypothetical" sources hinges on over a dozen 
spatial associations of unidentified EGRET sources at 
moderately low galactic latitudes with well-studied ra- 
dio SNRs. These include IC 443, 7 Cygni and W44, and 
many have neighbouring dense environments, which seem 
essential in order to provide sufficient gamma-ray lumi- 
nosity to exceed EGRET's sensitivity threshold. Confir- 
mation of these associations will probably only come with 
better gamma-ray angular resolution (e.g. at the soft 
gamma-ray band addressed by the INTEGRAL mission) 
and/or positive TeV detections. The remnants associated 
with several of the EGRET unidentified sources show an 
apparently low level of TeV emission, as determined by 
Whipple (Lessard et al. 1995), which may be marginally 
inconsistent with the spectra generated in the model of 
Drury, Aharonian and Volk (1994); Gaisser, Protheroe 
and Stanev (1996) use a slightly steeper accelerated pro- 
ton distribution to provide greater consistency between 
the EGRET and Whipple bands. De Jager and Mas- 
tichiadis (1996) suggest that there may be an intrinsic 
cutoff in the SNR-generated cosmic ray population that 
would spawn a complete absence of TeV emission, which 
is quite compatible with the data. 

All of the above models invoke simple power-law accel- 
erated particle populations. In this paper, we utilize the 
more sophisticated output of shock acceleration simula- 
tions (e.g. Jones and Ellison 1991) to address the issues 
of spectral curvature and the maximum energy of accel- 
eration in the context of SNR gamma-ray emission. We 
use output from the fully non-linear Monte Carlo simula- 
tions of Ellison and Reynolds (1991) and Ellison, Baring 
and Jones (1996) to describe the accelerated population 
in environments where it infiuences the dynamics of the 
SNR shell. Such non-linear effects are crucial to the mod- 
elling of gamma-ray emission in SNRs. Our results make 
clear predictions of what maximum energies of gamma- 
rays are expected, and further make more accurate de- 
terminations of the level of TeV emission relative to the 
MeV range in these sources. We find that the TeV/GeV 
fiux ratio in our non-linear regime is reduced by factors 
of a few below pure power-law proton scenarios, a result 
that is very important for correlated observations by IN- 
TEGRAL, Whipple and other gamma-ray instruments. 



2. FERMI ACCELERATION AT SHOCKS 



Diffusive shock acceleration is usually assumed to gener- 
ate power-law particle populations in astrophysical mod- 
els. This simple approximation omits the effect the ac- 
celerated particles themselves have on the hydrodynam- 
ics of their shocked environment. Considerations of such 
non-linear dynamics are essential to the gamma^ray SNR 
problem since the peak luminosity arises at the onset of 
the Sedov phase (Drury, Aharonian and Volk 1994) where 
the fraction of ram pressure going into cosmic rays and 
nonlinear effects are maximized. These non-linear effects, 
well-documented in the reviews of Drury (1983) and Jones 
and Ellison (1991), have a feedback on the acceleration 
mechanism and its efficiency. This becomes evident when 
it is observed that the "test-particle" power-law slope de- 
pends purely on the compression ratio of flow speeds on 
either side of the shock, and that this ratio is ultimately 
dependent on the hydrodynamics. 

Monte Carlo Proton Distribution 
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Figure 1; A typical proton distribution np{ep) , of arbi- 
trary normalization, resulting from Monte Carlo simula- 
tions of particle acceleration at SNR shocks. The transi- 
tion to relativistic energies produces a "bump" at around 
the kinetic energy Ep \ GeV, a result of assuming 
a momentum- dependent mean free path. The £'^np{£p) 
representation gives a good feel for the spectral curvature 
produced by the non-linear hydrodynamics above 1 GeV. 

Our Monte Carlo simulation of Fermi acceleration deter- 
mines the acceleration efficiency and maximum energy 
of acceleration. Although steady-state, this method can 
roughly mimic the net features of the time dependence of 
the SNR expansion and may describe the general spec- 
tral properties quite well just prior to the start of the 
Sedov phase. A typical proton distribution is shown in 
Fig. 1. Positive curvature appears in the distribution 
above 1 GeV due to the non-linear effects. The differ- 
ence in slopes about the proton rest mass energy at ~ 1 
GeV is just due to the varying dependence of momentum 
on energy as prescribed by special relativity. The pres- 
sure of the accelerated population acts to slow down the 
fast-moving flow upstream of the shock. 

High-momentum protons have longer mean-free paths A , 
a model assumption that is supported by observations 
made by plasma experiments at the Earth's bow shock 
(e.g. Ellison, Mobius and Paschmann 1990), and there- 
fore typically influence the flow on larger scalelengths. 



The maximum compression ratio is achieved at the largest 
scales, i.e. farthest from the downstream region (inside 
of the SNR shell; see Ellison, Baring and Jones 1996 for 
typical velocity profiles), implying an increased efficiency 
of accelerating higher-energy particles relative to those 
at lower energies. Upward spectral curvature ensues, 
and this has a siginficant influence on predictions of TeV 
fluxes in gamma-ray SNRs. Non-linear predictions of the 
spectral index differ significantly from the test-particle 
case (Ellison, Baring and Jones 1996); in Fig. 1, a 
proton momentum distribution from a strong shock would 
be flat (i.e. horizontal) above 1 GeV in the e^npi^Sp) rep- 
resentation. Deviations from a power-law behaviour can 
be up to a factor of around 3 over four decades in particle 
energy in the 1 GeV 10 TeV range. Such enhancements 
are extremely important for discussion of the potential 
observability of EGRET sources by Whipple, HEGRA 
and other ground-based TeV experiments. 



2.1. The maximum particle energy 



The maximum energy of Fermi-accelerated protons is de- 
termined by two considerations, namely the acceleration 
time vs. the age of the remnant, and the diffusive scale- 
length of the particles vs. the size of the ionized region. 
For test-particle shocks (i.e. linear ones that produce 
power- law protons), the acceleration timescale (e.g. see 
Forman and MorfiU 1979, Drury 1983 for derivations) for 
the Fermi mechanism at plane-parallel shocks up to en- 
ergy Et^y (in units of TeV), far above the injection en- 
ergy, can be written as 



31.7 



QB 



V r(r -I- 1) Et^y 



(1) 



*1000 



where r = Wups/wdown is the compression ratio, i?_5 
is the field in units of 10~^ Gauss, and r] is the ratio 
of the particle mean free path A to its gyroradius rg . 
Here Q is the charge number of the ion (e.g. Q = 1 
for protons and Q = 2 for alpha particles). Eq. (1) is 
a specialization of Eq. (4) of Ellison, Baring and Jones 
(1995), expressed in units appropriate to the SNR prob- 
lem at hand. Quite generally, r] > 1 and one infers 
values of 1 < 77 < 10 in interplanetary shocks in the in- 
ner heliosphere (Baring et al. 1997). Shock speeds itiooo 
( = Uups ) are in units of 1000 km/s. Implicit in such a 
formula is the assumption that the mean free path scales 
as X oc Vg oc p , which seems reasonable in the light 
of the aforementioned observations of plasmas in the he- 
liosphere. Although Eq. (1) is strictly confined to test- 
particle regimes, it serves also as an order-of-magnitude 
estimate for the full non-linear problem that is treated 
by the Monte Carlo simulation. Assume that the age of 
remnants is just their shell radii, typically in the vicinity 
of a few pc, divided by the expansion speed uiooo • This 
contention can be modified by the dynamical interaction 
of the supernova ejecta with the ISM, leading to larger 
age estimates (e.g. Chevalier 1982). For SN la, typically 
Wiooo 10 , while for the more-massive ejecta associated 
with SN II, the shock speed is slower: Miooo ^ 4 . Hence 
(bright) remnant ages are typically 1 — 10 X 10"'^'^ sec, 
implying age-limited acceleration terminating at around 
10* - 10^ GeV. 

The diffusion scale of the SNR medium is k/u for a dif- 
fusion coefficient of K = Xv/3 . Hence the diffusion scale 



(in parsecs) is comparable to 



rfpc = 3 X 10 
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For acceleration regions of the order of a tenth of the 
shell size (much smaller might be expected for remnants 
in the neighbourhood of dense neutral regions, which is 
frequently the case for remnants associated with uniden- 
tified EGRET sources), size- limited acceleration termi- 
nates also at around 10^ — 10^ GeV. The factor of a tenth 
represents an estimate of the effective escape lengthscale 
taking into account shell geometry and the remnant his- 
tory. Note that if the acceleration time in Eq. (1) is 
multiplied by the shock speed, the resulting length scale 
is of the order of Eq. (2), differing only by factors in- 
volving the compression ratio. This is not surprising 
since both estimates have their origin in particle dif- 
fusion. Clearly higher compression ratios will enhance 
the acceleration time, defining the general property that 
stronger shocks usually produce size-limited rather than 
time-limited maximum energies of acceleration. Another 
important property of these estimates is that Eqs. (1) and 
(2) depend on the charge of the species but not explic- 
itly on the mass. licnce electrons potentially have similar 
maximum energies to protons if cooling is absent, while 
helium and heavier nuclei with higher charge states can 
be accelerated to higher energies (as opposed to energy 
per nucleon). 



3. GAMMA-RAY EMISSION SPECTRA 



The effects of using the full non-linear Monte Carlo sim- 
ulation for predicting gamma-ray cniissivities from SNRs 
can be adequately demonstrated using just the tt*^ emis- 
sion component. Treatment of other radiation mecha- 
nisms (discussed below) is deferred to future work, in 
part because the hadronic contributions tend to provide 
the major portion of the gammarray radiation. The tt^ 
emissivity and spectrum were calculated much along the 
lines of the work of Dermer (1986): pp collisions produce 
tt" s, which subsequently decay to produce two photons. 
The proton component of the cosmic rays collides with 
nuclei in the cold ambient ISM; note that cosmic rays of 
higher mass number also contribute significantly to the 
TT emissivity. Low-energy protons, typically with kinetic 
energies below a few GeV, create pions via the A res- 
onance, following the model of Stecker (1970). As their 
energy increases, more resonances are sampled. High en- 
ergy protons, typically above 10 GeV, interact according 
to a radial scaling empirical formalism (e.g. Tan and Ng 
1983). Details of these formalisms, kinematics and nu- 
merical procedure, can be found in Baring and Stecker 
(1997, in preparation). Note that the gamma-ray spec- 
tral profiles for different cosmic ray proton kinetic ener- 
gies are symmetric about the energy rxiT^c? jl , due to the 
pion decay kinematics (e.g. Stecker 1970). 



Fig. 2 shows various representations of a tt'^ gainnia- 
ray spectrum resulting from a proton distribution gen- 
erated by the nou-lincar Monte Carlo simulation. In this 
case, protons are accelerated out to a few tens of TeV, 
a limit obtainable from Eq. (1) for a very young rem- 
nant with 77 ~ 10 and r ~ 11 (it must be remem- 
bered that non-linear shocks generate overall compres- 
sion ratios much larger than 4; see Jones and Ellison 
1991; Ellison, Baring and Jones 1996). The differential 



photon spectrum is largely visually uninteresting, so inte- 
gral and double-integral ( " I'Fi^ " ) spectra are plotted to 
highlight the spectral properties. One prominent feature 
is that spectra flatter than result. The spread in- 

duced by the kinematic phase space for pion production 
and decay tends to smear out the underlying curvature of 
the proton population. Yet, for comparison, the integral 
spectrum resulting from a power-law proton population 
with a high-energy cutoff is depicted in Fig. 2, indicat- 
ing that it overestimates the TeV/EGRET flux ratio by 
a factor of a few. This clearly yields an indication of 
the improvement of the non-linear calculation over the 
standard test-particle p^^ infinite power-law case con- 
sidered by Drury, Aharoiiiaii and Volk (1994, and also 
approximately by Gaisser, Protlieroe and Stanev 1996). 
Therefore it follows that these test-particle implementa- 
tions of shock acceleration theory can significantly over- 
predict the TeV/EGRET flux ratio. Consequently, the 
incorporation of non-linear dynamical considerations can 
potentially relax observational constraints imposed by the 
Whipple upper limits to IC443, W28 and 7-Cygni as pre- 
sented by Lessard et al. (1995). Note that the spectral 
structure below 1 GeV in the proton distribution in Fig. 1 
is immaterial since such kinetic energies are below the 
production threshold. 



Gamma— Ray Flux 
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Figure 2: The -^-ray flux at earth, for a young SNR about 
a kpc distant, integrated over a shock acceleration pro- 
duced proton distribution (long dashed lines). The solid 
curve represents the integral spectrum for this distribu- 
tion, and the dash-dot curve the integral spectrum from a 
pure power-law ( p~^ ) proton population ( exponentially 
cutoff at around 1 Te V) with index given by the low en- 
ergy portion of the non-linear model proton distribution. 
The ei^n^ representation clearly illustrates the effects of 
spectral curvature and the cutoff. 

DeJager and Mastichiadis (1996) suggested that con- 
straints to the TeV/EGRET flux ratio can be relieved by 
having the photon spectrum turn over below the Whipple 
sensitivity threshold. The maximum energies of the pro- 
tons at ^ 10 — 30 TeV propagates down into the few TeV 
range in the gamma-ray spectrum due to the inelasticity 
of the pp ppir'^ — > PPll chain. Hence for the illus- 
trated case, HEGRA would be unlikely to see an EGRET 
source, but Whipple potentially could. Turnovers at such 
low energies are predicted neither in the work of Gaisser, 
Protheroe and Stanev (1996), nor by Drury, Aharonian 
and Volk (1994), mainly because they focus on cooling- 
limited rather than age/size-limited acceleration. Their 
inverse Compton cooling rates off the microwave and local 



IR backgrounds turn out to be low enough that the freely 
expanding shell purely defines the maximum proton en- 
ergy. The non-linear shock acceleration solution treated 
here provides natural cutofTs in the TeV range that can 
also yield compatibility with the Whipple upper limits for 
the appropriate choice of remnant shock parameters. By 
inspection of Eq. (2) , slightly lower shock speeds can sup- 
press emission in the TeV range, which may well be nec- 
essary to explain the four orders of magnitude difference 
in integral flux inferred for 7-Cygm (see Lessard et al. 
1995) between 100 MeV and 500 GeV from the EGRET 
source flux and the Whipple upper limits. The implica- 
tion of such low maximum energies for cosmic ray protons 
is naturally that older SNRs arc required to produce the 
bulk of cosmic rays out to the lO^^ey "knee," so that 
gammarray SNRs are a gamma-ray bright minority of the 
cosmic ray-producing remnant population. 



4. DISCUSSION 



Generalizing beyond 7r*^ emission generated via pp colli- 
sions is clearly important. These interactions also spawn 
charged pions at comparable rates, and these yield sec- 
ondary electrons and positrons and consequently inverse 
Compton emission as discussed in Gaisscr, Protheroe and 
Stancv (1996). While the tt'^ decay spectrum traces the 
proton distribution, the inverse Compton spectrum is po- 
tentially much flatter due to the kinematics. Hence, in 
principal, the inverse Compton mechanism can be much 
more efficient at higher energies, and perhaps dominate 
the tt" decay products; such scenarios arc discussed in 
the context of the difi'usc galactic gamma-ray background 
by Hunter et al. (1996). Primary shock-accelerated elec- 
trons are in principal possible, though the theoretical ef- 
ficiency of their generation is very poorly understood. 
Gaisser, Protheroe and Stanev (1996) argue that the 
EGRET/Whipple spectral constraints for IC 443 imply 
that inverse Compton emission is relatively unimportant 
and also that the primary electron population is of low 
abundance. Within the confines of non-linear shock ac- 
celeration theory, unlike the power-law choices made (i.e. 
linear shock models) by Gaisser, Protheroe and Stancv 
(1996) and De Jager and Mastichiadis (1996), the accel- 
erated electron spectrum does not trace the proton dis- 
tribution at all energies. Below a few GeV, the electron 
diffusion lengths are always less than those of protons of 
comparable energy, so that electrons effectively experi- 
ence a weaker portion of the shock structure than do the 
protons (e.g. see Ellison and Reynolds 1991 for a dis- 
cussion). The net effect is that the electron distributions 
arc steeper below a few GeV and the population level 
above a few GeV much lower than is predicted in the 
linear shock models. Therein lies another crucial aspect 
of the inclusion of the non-linear effects of shock accel- 
eration theory: the linear approaches cited above seri- 
ously overestimate the efficiency of electron acceleration, 
and hence the contribution of the inverse Compton and 
bremsstrahlung processes to the SNR luminosity. Such 
issues will be dealt with in a future paper (Baring et al. 
1997, in preparation). 

The issue of the prominence or otherwise of an inverse 
Compton component may well be elucidated by coordi- 
nated observations at the TeV range and by INTEGRAL 
at MeV energies: the 7r° decay channel produces much 
lower emissivity below m,rc'^/2 ~ 67 MeV due to decay 
kinematics, so that an inverse Compton spectrum can be 
much more visible in the soft gamma-ray regime. An 



asset of INTEGRAL will be its angulax-resolution capar 
bility (e.g. see Winkler 1996). EGRET has difficulty 
resolving on scales of the order of the angular diameter 
of nearby remnants; hence it is not clear whether any de- 
tected gamma-ray emission is associated with whole rem- 
nants, portions or "hot spots" (or neither). Improved an- 
gular resolution at MeV and TeV energies will clearly pin 
down or disprove these associations. Note also that, apart 
from protons, heavier nuclei make up the cosmic ray pop- 
ulation and the ambient ISM, and these also contribute 
significantly to the gamma-ray and electronic products; 
their inclusion in our gamma-ray emissivity calculations is 
deferred to future work. In conclusion, coupled TeV/sub- 
GeV observations of supernova remnants, in which the 
INTEGRAL can play a big role, will discriminate between 
various gamma-ray emission models of these sources, and 
enhance our understanding of the cosmic ray production. 
The non-linear effects described here are an essential in- 
gredient for any model that invokes Fermi acceleration 
at SNRs, particularly since their prediction of enhanced 
TeV/GeV flux ratios can tighten model constraints im- 
posed by current TeV upper limits and possible future 
detections. 
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